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1. Introduction 

Chiral perturbation theory for Np = 3 light flavours at next-to-leading order was 
shown by Gasser and Leutwylerjl[ to involve ten independent coupling constants 
L[(/x) (referred to, sometimes, as LEC's). In ref. all ten couplings were deter- 
mined from low-energy data except two, L 4 and L§. According to the OZI rule or, 
alternatively, large N c considerations these couplings are expected to be suppressed 
relative to the other ones0. Physically, L 4 and L 6 control how observables of the 
pion, like M n or chiral SU{2) order parameters (F n )su{2), (uu)su{2) vary if the mass 
of the strange quark mass m s varies. If OZI suppression holds these quantities are 
expected to be essentially insensitive to variations of m s and the value of the chiral 
order parameters F^, (uu) in the 577(2) chiral limit or in the SU(3) chiral limit (i.e. 
m s = 0) ought to be very nearly the same. In nature, of course, the values of the 
light quark masses are fixed but they can be made to vary in unquenched lattice 
simulations of QCD. 

There are reasons to suspect that the couplings L 4 and L 6 are, in fact, far 
from being suppressed. One reason is suggested by two recent unquenched lattice 
simulations which have investigated the phase structure of Q CD-like theories with 
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Np equal- mass flavours when the value of Np is varied. It was found in ref.0 
that for Np > 6 a phase with no confinement and no spontaneous chiral symmetry 
breaking prevails, at any value of the QCD coupling constant. In ref.[§] a very strong 
decrease of the chiral order parameters was observed upon varying Np from Np = 2 
to Np = 4. The values of the coupling constants Li encode informations on the 
physics of the massive states in QCDQ]. In particular, L4, Lq are linked to the scalar 
resonances. It is notorious that the OZI or large N c rules seem to fail in this sector. 
In fact, is not clear at present exactly which of the scalar resonances form the lowest 
lying nonet (see e.g. the reviews in the pdg0). In other words, the fact that L 4 , L 6 
may be unsuppressed, the failure of the OZI rule in the scalar sector are connected 
and these features may be related in an interesting way to the phase structure of 
QCD, in particular to the fact that a phase transition could occur for a value of Np 
not exceedingly larger than Np = 3. An interpretation of these features in terms of 
a paramagnetic effect of the quark loops is discussed in ref . || 

These topics have started to be investigated in a previous paper 0. The couplings 
L 4 , L 6 were related to the scalar resonances via the scalar form-factors of the pion 
and the Kaon. These form-factors can be reconstructed directly from experimental 
data on S-wave scattering, modulo a few plausible hypothesis^]. The coupling L4 is 
related to the derivative of the strange form-factor of the pion at the origin and, for 
Lq, a chiral sum rule can be derived in the terms of the correlator Il 6 of the strange 
scalar current ss and the non-strange current uu + dd. In the present paper, we 
investigate the 0{m s ) corrections to the results of ref.0, which could be sizable. For 
this purpose, one needs to use CHPT at order p 6 . Renormalizability of the theory 
at this order was recently proved||, and the set of independent chiral lagrangian 
terms classified [[H]]. In particular, we will use the expansions of M n and F n which 
were derived in refs.|ll[] and [0 and we have computed the 0(p 6 ) expansion of the 
correlator Il 6 . The result of this effort, at first sight, will turn out to be mitigated: 
as new coupling constants appear in the formulas, there is a loss in predictivity over 
the leading order calculation. As discussed in sec. 2, conclusions may nevertheless 
be drawn provided one makes some assumptions on the convergence of the chiral 
expansion. 

In order to gain further insight, we have also investigated a different, more model 
dependent approach in sec. 3. Following the idea of ref.|| one starts from a lagrangian 
for the scalar resonances and the LEC's are generated by integrating out the reso- 
nances. If we attempt to determine some 0(p 6 ) LEC's in this way, from the most 
general lagrangian, we find that there are many undetermined resonance parameters 
which get involved, like three resonance couplings, couplings of one resonance to two 
scalar sources etc... In sec. 3, we discuss the more specific dynamics generated by the 
linear sigma- model [13]. This model is still attracting interest WM [T5[ in connection 



with the phenomenon of disordered chiral condensate [^] which could be formed in 
heavy ion collisions. In the symmetry breaking sector, we will consider a lagrangian 
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more general than previously done, still compatible with the criterion of softness (re- 
stricting ourselves to first order symmetry breaking). The model can accommodate 
exactly, in principle, a given set of scalar nonet masses and one obtains the 0(p A ) 
and 0(p 6 ) LEC's as definite predictions. The question at this point is which scalar 
mesons are to be included in the nonet? In particular, does one have to include a 
light a and a light k meson? It has long been known that an S-matrix pole can be 



identified in low-energy tttt scattering [18] with a very large imaginary part. A similar 



structure was argued to exist also in ttK scattering |]19| |20| , but it is unclear whether 
such poles should be interpreted as physical scalar resonances. For our purposes, we 
will leave this question open and consider two different possibilities for the scalar 
nonet. As we will see, in the framework of the sigma- model, these will correspond 
to rather different behaviour of the chiral vacuum. 

2. L4, Lq beyond 0(p A ) 

Our estimates of L 4 and L 6 are based, essentially, on a method to extract the scalar 
form-factors of the pion and the kaon using experimental tttt — KK scattering data, 
proposed in ref.||. The form-factors get determined up to a normalization factor, 
for which one uses CHPT. In order to determine L 4 one equals the experimental 
determination of G' n (0), the derivative of the strange scalar form-factor of the pion 
and the chiral expansion of this quantity. In ref. |7]] this matching was performed at 
leading order in the chiral expansion, i.e. G ; n (0) was expanded up to 0(p A ) and for 
the normalization condition we used CHPT at order p 2 . In a similar way, for L 6 we 
use the scalar form factors to express the spectral function of the correlator ^(s) 
(see (|2.13| ) below) and L 6 is obtained by matching I1 6 (0) evaluated from experimental 



data and its chiral expansion. Below, we discuss the chiral corrections to these results, 
which involve two parts a) one must use the chiral expansions of G' n (0) and II 6 (0) up 
to 0(p 6 ) and b) we must use CHPT at 0(p A ) instead of 0(p 2 ) in the normalization 
conditions. 

Let us designate by F n (s) and G n (s) the non-strange and strange scalar form 
factors of the pion, and by Fk{s) and Gk{s) the analogous ones for the Kaon. We 
will use the following normalizations, 

1 /3/nL-7„. 1 Xji_o_o\ n — 1 . n\ 



F n = jf o ]/%(0\uu + dd\7i°7r°) G n = — y-(0|ss|7rV> 

F K = 4-V2(0|«u + dd\K + K-) G K = -^V2(0\ss\K + K-) . (2.1) 

The method of ref. || consists in solving numerically a set of Muskhelishvili-Omnes 
coupled-channel equations using experimentally determined tttt — KK S-wave scat- 
tering T-matrix elements as input (in practice we made use of the parametrisations of 



ref.[]2T| and ref.p2|). Strictly speaking, these equations hold under the assumption of 
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exact two-channel unitarity up to s = oo while in practice, two-channel unitarity is a 
good approximation up to the 7777 threshold. The resulting form factors are expected 
to be reliable in a finite energy range, which we will assume to extend up to 1 GeV. 
One also chooses appropriate boundary conditions for the T-matrix at s = 00 which 
insure existence of a solution with a minimal number of free parameters: a solution 
vector u\(s), 7/2 (s) will be determined in the entire energy range once the values at 
one point, say at s = 0, are given. Extension to more channels and the stability of 
this scheme were discussed in ref. 0. In practice, one first constructs numerically 
two independent solutions of the equation set: it»(s) and Vi(s), i — 1,2, normalized 
at the origin such that iti(0) = 1, 1*2(0) = and t>i(0) = 0, ^(O) = 1. The form 
factors are then given as 

(2 - 2) 

and similarly for the strange form factors. In refs. || [[/J the normalization at the origin 
was taken from CHPT at 0(p 2 ). Here, we wish to investigate the chiral corrections 
to these results, so we need to go to the next chiral order. 

2.1 Evolution of F T 

In the SU(2) chiral limit, the strange form factor of the pion vanishes at the origin 
Cr,r(0) = 0, the LEC L 4 is related to the first derivative at the origin of G n . It is 
convenient to consider the proportional quantity 

d e ; p = sj^m s B Q G'M = ^m s B Q G K ^)v'M . (2.3) 

Here, v[(0) can be determined from experiment by the procedure outlined above, 
and one finds, 

v[(0) ~ 0.27 GeV" 2 , (2.4) 

using the T-matrix parametrisation of ref . (using that of ref . ]21| one would find 
v[(0) ~ 0.31 GeV" 2 ). We next need to express m s B and G K {0) in eq.(fD|. At 
leading chiral order, one has m s B Q = m? K — m 2 /2 and G K {0) = y/2. Including 
0(m s ) corrections to this result, one obtains d e p P in the form, 



TUk ( „ ~ 1 



m 2 K - mi/2 + ^-8S 8 - 16S 6 + — (L„ + 1] 



«i(0), (2-5) 



with 



2 
in 

S 8 = -2L r M + L r M, S, = -2L r M + L r M, L v = log . (2.6) 

/i 

Corrections of order rri^m 2 K or m\ have been neglected here and CHPT at 0(p A ) 
has been used to obtain Gk{0) and to express m s Bo in terms of the physical Kaon 
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mass. We observe that low-energy constants appear now but it its consistent to use 
their values obtained at 0(p 4 ), i.e. for L 5 , L 8 the values given in ref.[[J and for L 4 , 
Lq those obtained in ref.jTj. We can then match d e p P with its expression as a chiral 
expansion. The latter is obtained from the relation || 

m s dF n 

dF = ~77~~, — 2.7 

F n dm s 

which holds in the chiral SU{2) limit, m u = = 0, and one can use the chiral 
expansion of F n which has been determined up to 0(p 6 ) [fLl]|[]12" |. One obtains 

4=^ + ^+0^), (2.8) 



with 



and 



d< " = 8L ' ~ 3^ {Lk + 1)1 iA ' = log! ^ £ • (2 ' 9) 



7T" 



vr 2 V 768vr 2 9 1 9 2 9 3 9 V V ; 

1 /26 35 5 3 3 0.0022X 

+ ^ U Ll + 54 La + 9 L4 + 4 5 ~ 3Le ~ 2 8 " • 

In this formula, the finite contribution from the so-called sunset diagrams has been 
evaluated numerically. Equating (|2.3|) and (|2.8|) at the leading, linear order in m s , 
gives the 0(p 4 ) determination of L 4 . From the numerical results of refs.||0 one 
obtains 

Ll(m v ) ~ 0.6 10" 3 [order p 4 } . (2.11) 

If we include the corrections quadratic in m s now, a new 0(p 6 ) coupling constant 
appears, labelled Ci 6 in ref.|TlJ, so strictly speaking, we have one equation and two 
unknowns. A similar situation will prevail in the determination of L§ to be discussed 
below. We expect, however, the chiral expansion to be meaningful and the 0(p 6 ) 
part of the expansion to be smaller than the 0(p 4 ) one. If we assume a given ratio for 
these two parts we can determine L 4 and C\§ separately. Once these two couplings 
are known we can, furthermore, using the CHPT expressions of ref . [fLl]! Ul2| , determine 
how F n varies in going from an SU(2) chiral limit to an SU(3) one. Some results are 
collected in table 1, where we have varied the 0(p 6 ) over 0{p 4 ) ratio between 10% 
and 100%. 

This exercise seems to indicate that the determination of L 4 is rather stable and 
rather close to its determination at leading order, and provides an estimate for the 
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0(p 6 )/0(p 4 ) 


10 3 L 4 


io 5 c 16 


F/F - 1 


0.10 


0.58 [0.31] 


0.36 [0.24] 


0.18+0.005 


0.30 


0.53 [0.26] 


0.39 [0.26] 


0.17+0.011 


0.50 


0.49 [0.22] 


0.41 [0.28] 


0.16+0.016 


1 


0.42 [0.15] 


0.44 [0.32] 


0.14+0.024 



Table 1: Assuming a given ratio of the 0(p 6 ) to the 0{p A ) contributions in the chiral 
expansion of cLf, the values of the LEC's L\{n) and C[ 6 (/i) are computed for fi = m v 
(in brackets, = m p ). We also display the successive contributions to the ratio of F = 
(F 7T )su(2) over F = (F n ) su{3) 

size of the 0(p 6 ) coupling constant C[ 6 . A priori, however, one cannot completely 
exclude a different solution, with L 4 (m J? ) = 0, for instance. In this case, all the 
OZI violation would be concentrated in the 0(p 6 ) parameter, which seems somewhat 
unplausible. Taking into account the experimental uncertainties in the T-matrix, the 
uncertainty from the energy region above 1 GeV (which were estimated in ref.0) and 
the discussion above of the 0(m s ) corrections, we find that the error on L4 should 
be of the order of 30%, i.e., 

£ 4 (m„) = (0.6 ± 0.2) 10~ 3 . (2.12) 



2.2 Evolution of (uu) 

Let us now discuss the coupling constant L 6 . Using the same scalar form factors, it 
was remarked in ref. that one can derive an estimate of Lq. The idea is to consider 
the two-point correlation function, 

MP 2 ) = -4 / d 4 xe ipx (0\T(uu(x) + dd(x))ss(0)\0) ^ . (2.13) 

In an energy domain extending roughly up to the f]f] threshold, the spectral function 
is given in terms of the scalar form factors of the pion and the Kaon, 



+ \/ ^F K (s)G* K (s)9(s - Am 2 K ) ■ (2-14) 

Essentially, a usual assumption is made there, that the 4n or Qn contributions are 
negligibly small in this energy region (see, e.g. ||21||). In the 577(2) chiral limit, m u = 
md = 0, a super- convergence relation holds, 



ImU 6 (s)ds = . (2.15) 
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Similar chiral sum rules are known to be saturated to a fairly good approximation in 
terms of a few low-lying resonances (see, e.g.||23||). In the energy region below 1 GeV, 
inserting the form factors constructed as discussed above, the contribution to the 
integral is found to be positive, and dominated by the / (980) resonance. A plausible 
assumption, then, is that the higher energy contribution to the sum rule ( |2.15| ) will 
be saturated, at least approximately, by the next prominent scalar resonance, the 
/o(1500), which we expect to make a negative contribution. We observe that the 
signs of the contributions of the /o(980) and the /o(1500) conforms with the recent 
assignments of Minkowski and Ochs[]24| as (essentially) SUf(3) singlet and octet 
respectively. Finally, one can calculate the value of the correlation function at zero 
energy, n 6 (0), making use of the scalar form factors below 1 GeV, and the super- 



convergence relation ( |2.15| ) as 



n 6 (0) ~ - ( f 1 /mIl6(g) ds-- f 1 7mn 6 (s) ds) , (2.16) 

7T \J0 S S JO J 

with so — 1.5 2 GeV 2 . This gives us an "experimental" value of IT^O), the difference 
with ref. [/J is that we now wish to take into account chiral corrections in the normal- 
ization of the form factors. Below the KK threshold, we can express the spectral 
function in terms of the basis solutions to the Muskhelishvili-Omnes equations Ui(s), 
Vi(s) as, 



ImU 6 (s) = F n (0)G K (0) Ul (s)vl(s) 

+F K (0)G K (0)v 1 (s)v* 1 (s)] . (2.17) 

Including 0(m s ) chiral corrections in the normalizing factors, we now have 



F n (0)G K (0) = 2y/3 1 + 



m 



K 



F K (0)G K (0) = 21 + 



m K 

F2 




2 r 7 



32S 8 -80S 6 + —L v + —_), (2.18) 



where ^8 = — 2L§ + L5, ^6 = — 2Lg + L\ and we may use in this part the values of 
the LEC's as determined at 0(p 4 ). 

We then use the experimental result for n^O) in conjunction with its chiral ex- 
pansion. The computation of Hq(s) to 0(p 6 ) is exposed in the appendix. Concerning 
I1 6 (0), the chiral expansion goes as follows 

n 6 (o) = nw + ^n( 6 ), (2.19) 



F 2 



with 

1 /22 



n(4) = 64 ^-i6^(y (L - + 1) + 9 L43 ) ' (2 - 20) 
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and 

L x = log^/, L 43 = log-. 2.21 
\i A 3 

This expression generates the 0(p 4 ) determination^] of L 6 , 

Ll(m v ) ~ 0.5 1(T 3 [order p 4 ]. (2.22) 
The expression of the 0(p 6 ) part in ( |2.19| ) is, 



2 



J_T2_J_ 

72vr 4 K 967r 4 ^ 43 ^ 288vr 4 " 43 
L K /64 70 16 55 



n( 6 ) = 256(C 2 r n + 3C") T Li rL^L 43 + -L 

+ 



7T 2 V 3 9 9 25927T 2 ; 

W 16 9 + 16 9 + 16 S 47 ^ (223) 
'IF \~3 S6 + Y Ss + T 7 ~ 5184^ J (2 ' 23) 

1 /140 175 8 19 \ 

V I 9 6 + "^T 8 + 27 7 ~ 5184^ J ' 



with 



S 7 = 3L 7 + Lg. (2.24) 

In these formulas, we have set m u = ma = 0. These expressions exhibit the explicit 
m s dependence (which will prove usefull below as we plan to integrate over m s ). For 
numerical application, we may use m s B = m 2 K — m 2 /2 in n( 6 ) while in we 
have to use the 0(p 4 ) expansions of m 2 K and m^jl]]. Equating the chiral expansion of 
ne(0) with its sum rule evaluation gives us a linear relation involving the 0(p A ) LEC 
L 6 and the 0{p & ) combination C 20 + 3C21. As before, we must make a hypothesis 
concerning the convergence of the chiral expansion if we want to separately estimate 
these two contributions. 

Another interesting information contained in the correlation function Tl§ concerns 
the variation of the quark condensate in the chiral SU(2) limit as a function of the 
strange quark mass: 



d(uu) B t 



IT 6 (0) . (2.25) 



dm s 2 

Using the chiral expansion of 116(0) we can integrate this equation in the variable 
m s from its physical value down to m s = 0. In this way we obtain an estimate for 
how the condensate varies from a chiral 577(2) limit to a chiral SU(3) limit. This 
variation is obtained as an expansion in powers of the physical strange quark mass 
(or, alternatively, in powers of the Kaon mass), 

(uu) sm = (uu) sm (l + + ^R^) , (2.26) 



with 

^ 4) - 32 ^-144^ L --72^ L - (2 ' 27) 
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and 



7T 2 V 3 9 18 10368tt 2 
L 43 /4 „ 4 „ 4 „ 5 

vr T V3~ u ' 9"' ' 9"° 51.S17T 
H — n ( "5*6 — —5*7 + —5*8 



^ U* + o*+o*-»i7I5) (2 - 28) 



vr 2 V 9 ° 27 ' 54 
Using the relation 

-2(n M ) st/(2) = ^-(m^ 2 ) m=o , m =I(m M + m ,), (2.29) 

one sees that the expansion of the condensate can be rederived from the expansion 
of the product Tn^F 2 in powers of m s . We have verified, using the 0(p 6 ) expansion 
P@ of m n , F n that the formulas ( gg§ , (F37D , ( ggg ) are exactly reproduced. 

Let us now discuss some numerical results. At first, it is instructive to compare 
the value obtained for n 6 (0) using the normalization of the form factors at 0(p 2 ) 
and the normalization at 0(p A ), 

n 6 (0) ~ 0.022 [0(p 2 ) norm.], n 6 (0) ~ 0.043 [0(p 4 ) norm.]. (2.30) 

We see that including 0(m s ) corrections in the normalizations brings a rather sub- 
stantial change in the result by approximately a factor of two. This can be traced 
to the large numbers appearing in front of the combination Sq = — 2L 6 + L 4 in the 
normalization factors ( p. 18 ). Evidently, the difference — 2L 6 + L 4 is very sensitive 



to a small variation of L 6 or L 4 . Taking L 6 slightly smaller than the central value 
(|2.22 ) would decrease significantly the modification in rLj(0) without essentially al- 



tering the other results. Upon considering the chiral expansion of ILj(0) now, this 
result suggests that the 0(p 6 ) contribution could be substantial and it is likely to be 
positive. We have collected some results in table 2 assuming that the ratio of the 
0(p 6 ) to the 0(p 4 ) contributions is positive and ranges from 10% to 100%. In this 
range, we find that the rate of convergence for the ratio of quark condensates seems 
reasonable. 

Both and (uu) show a tendency towards chiral restoration in going from 
chiral 577(2) to chiral SU (3). This tendency seems much stronger for (uu). This is in 
agreement with the arguments of ref.|| based on the spectrum of the Dirac operator. 
One must however bear in mind that its dimensionality is different from F n : if we 
had considered (uu) 1 ^ 3 , we would have found a smaller variation. This behaviour 
could perhaps suggest the possibility, at larger N F , of a phase with vanishing quark 



condensate, yet with chiral symmetry still spontaneously broken f25|]. Kogan et al.p6| 



have discussed how this could result from a discrete Z2 axial subgroup remaining 
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0(p 6 )/0(p 4 ) 


10 3 L 6 


10 5 (C 20 + 3C 21 ) 


(uu)su(2)/(uu)su(3) ~ 1 


0.10 


0.73 [0.56] 


0.29 [-0.02] 


0.85 + 0.14 


0.30 


0.64 [0.47] 


0.37 [0.06] 


0.76 + 0.18 


0.50 


0.57 [0.41] 


0.43 [0.12] 


0.69 + 0.21 


1 


0.46 [0.30] 


0.52 [0.21] 


0.58 + 0.26 



Table 2: Assuming a given ratio of the 0(p 6 ) to the 0{p A ) contributions in the chiral 
expansion of lle(0) the values of the LEC's Lg(m„) and the combination C^oOu) + 36*21 (/■*) 
are computed for p = (in brackets, [i = m p ). We also display the successive terms in 
the expansion of (uu) su(2)/ '(uu) S u(3) ~ 1- 

unbroken, but argue against this possibility in QCD. Concerning the error in this 
evaluation of L 6 , finally, it is expected to be somewhat larger than the error on L 4 
because of a greater sensitivity to the energy region above 1 GeV. The discussion 
above also suggests that one should have L e (m v ) < L^ijn^) otherwise the 0(m s ) 
corrections could become out of control. Keeping this mind, we find that L 6 should 
lie in the following range, 

L & {m v ) = (0.5 ±0.3) 10" 3 . (2.31) 



3. The LEC's in a linear sigma- model 

3.1 General scalar meson lagrangian 

Here, we adopt a simple resonance saturation point of view for estimating the low- 
energy coupling constants, which was discussed in detail in ref. J4|. It consists in 
making a tree level calculation starting from a lagrangian for the resonances. We 
will consider the scalar resonances here. It is most convenient to start from a repre- 
sentation of the resonances in which they transform under a non-linear representation 



of the chiral group [^]. A detailed discussion can be found in ref.Q, we adopt es- 
sentially the same notations here. If one is interested in the 0(p A ) LEC's then one 
needs consider only those terms which are quadratic in the resonance fields, 

C ss = i + l(V»SV»S) - l -M%Sl - \Mi(S 2 ), (3.1) 

or containing one resonance field and one scalar source, or one scalar field and two 
chiral fields w M , i.e. 

£s x = c d (Su^) + CdSoiu^Ufj) + c m (Sx+) + c m S (x+), (3.2) 

where So is the singlet scalar, with chiral limit mass Mq, S is a traceless matrix 
encoding the scalar octet, which has a common mass Mg in the chiral limit. We 
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use the same conventions otherwise as ref.jp. This lagrangian yields the resonance 
saturation estimates for the following low-energy constants, 





+ 




^5 = 




3M| 


M 2 ' 


M|' 


c 2 




c 2 




c 2 


+ 


^8 = 




6M| 


2M 2 ' 


2M| 



(3.3) 



In addition, the scalar resonances make contributions to Li, L 



3, 

Ll = ~6M| + 2M 2 ' L3= 2M|' (3 ' 4) 



which receive other important contributions from the vector and axial-vector res- 
onance sector. These tree-level estimates yield coupling constants which are scale 
independent 1 : one usually assumes that they should represent meaningful estimates 
of the Ll(fi) for values of [i for which chiral logarithms are numerically small, that is 
of the order of fi = 0.5 to /i = 1 GeV. In ref.[|J the OZI rule was assumed to apply, 
implying L4 = Lq = 0. The formulas (|3.3|) provide one relation between the values 
of L 5 , L s and the experimental value of the decay width — > in]. Values of the 
couplings Cd and c m were extracted, 

\c d \ ~ 32 MeV \c m \ ~ 42 MeV . (3.5) 

Below, we will investigate the more detailed predictions that one can make if one 
assumes the validity of a renormalizable, sigma-model type interaction lagrangian. 
Before we discuss this, let us consider the generalisation of the considerations above 
to LEC's of chiral order 0(p 6 ). We will illustrate the complications which arise at 
this order by considering only the lagrangian terms which are cubic in the scalar 
source x+- These terms involve the three coupling constants [[ToJ C\g, C20, C21, 



= ... + C 19 <4> + C 20 (xl)( X+ ) + C 21 <X+> 3 + ... (3.6) 

Firstly, one needs to consider interaction terms which are cubic in the resonance 
fields, 

Csss = aS 3 +bS (S 2 ) + c(S 3 ), (3.7) 
then terms which are quadratic in the resonance fields and linear in the scalar source, 

£SS X =Zm (S ){X+)+ fm S (SX+) 

+ g° m (S 2 X+ )+h m S 2 (x + ), (3.8) 



x If one starts from the renormalizable sigma-model it is possible, in principle, to compute the 
one-loop effective action and then generate the 0(p 4 ) LEC's with the correct scale dependence. 
This was done for the SU(2) sigma-model in rcf.pq 
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and, finally, terms linear in the scalar field and quadratic in the scalar source, 



£s X x = a ' s o(x+) 2 + b'S (xl) + c'(Sx+)(x+) 

+d'{Sxl). ( 3 - 9 ) 

This already brings in a large number of couplings and it would not be possi- 
ble to make definite predictions without additional assumptions. The terms in 
eqs.(|3~T| ) ,( |3.8| ) and ( |3.9| ) generate chiral 0{jP) coupling constants as well as chiral 
symmetry breaking for the scalar meson masses. It is convenient, both for generat- 
ing the 0(p 4 ) and 0(p 6 ) LEC's and for expressing the scalar meson masses to make 
the following field redefinitions, 



The redefined fields So and 5* have, by construction, no trilinear couplings and no 
minimal one-resonance to one-source coupling. The couplings of two resonances to 
one scalar source, eq. (|3T8|) , get redefined to 



_o 3c m c , _,° 3c m a 

9m 9 m ' m m J^2 ' 

o c m 6 c m c o 2c m b , . 

e m =e m - j^, J m =j m +^ . (d.ll) 

The scalar meson masses, to linear order in the quark masses, can be straightfor- 
wardly expressed in terms of these redefined parameters. The masses of the isospin 
1 = 1 and 1=1/2 mesons read, 

M l = M l ~ 4e m (2m^ + ml) - Ag m ml 



M 2 Ko = M 2 ao - 4g m (m 2 K ~ O , (3.12) 



while the isospin zero meson masses diagonalise the singlet-octet mass matrix which 
gets parametrised as follows 

Mn = Ml - Ah m (2m 2 K + ml) 

M12 

M 22 = Ml - y g m (m 2 K - mi) . (3.13) 

In deriving these expressions for the masses, we have ignored possible terms in the 
scalar lagrangian of the form d fl Sod' 1 So(x+), (V fl SV ll Sx+) , anticipating on the fact 
that such terms will not appear in the model to be discussed below. 
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Concerning the LEC's Cxg, C 2 o, C 2 i, the following expressions are obtained, 
F ° Cl9 -Mi{ 9m+ Mi) +d Mi 

■p-2(~i C m ( ^ o _i_ CC m\ . C m C m f ° b C m \ 



F - 2 (C 20 + 3C 21 ) - (<? ra +^|J + (/i° m +^|) 

3.2 Parameters of the linear sigma-model 

In the SU(3) linear sigma-model, one first encodes the scalar nonet and the pseudo- 
scalar nonet into a complex 3x3 matrix E, which, under a chiral transformation 
transforms as, 

£ - g R T.gl , (3.15) 

and the lagrangian is assumed to be renormalizable. In the chiral limit, the most 
general renormalizable lagrangian, invariant under the chiral group (except for the 
{/a(1) subgroup) contains four parameters, 

C = \(d^d^) - ^ 2 <££t> _ A((££t) 2 > 

-A'^E 1 ") 2 + /3 (det(S) + det(S f )) . (3.16) 



This lagrangian was reconsidered recently by several authors [ 15| [29 ] [ 3C ] [ 31 ] and we 
use the same notations as in ref. |5T|. Suitable choices of the parameters /i 2 , A, A', f3 
ensure that spontaneous chiral symmetry breaking occurs, i.e. that the potential is 
minimised for a non-zero value of the vacuum expectation value of the £ matrix, 

£ = vl d , (3.17) 

(which we take to be diagonal, assuming that no spontaneous breaking of flavour 
symmetry occurs). Values of v which correspond to extremums of the potential are 
solutions of the cubic equation, 

v ((4A + 12A> 2 - 2pv + /i 2 ) = 0, (3.18) 

which has real solutions v ^ provided the parameters satisfy 

(3 2 > ^ 2 (4A + 12A') . (3.19) 
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Next, expanding around the minimum we find that the pseudoscalars are massless, 
except for the singlet, r] which mass is proportional to j3, 



M 2 m = 6(3v . (3.20) 

We can also read off the expressions for the singlet and octet scalar meson masses in 
the chiral limit, M and M 8 , 

Ml = 8Xv 2 + 4f3v, Ml = 8(A + 3A> 2 - 2(3v . (3.21) 

Stability of the vacuum requires that these squared masses be positive. Using 
eqs. ( |3.18|) ( |3.20|) and ( |3.21| ) we can trade the original four parameters of the lagrangian 



( |3.16| ) for the more physically relevant ones v, M vo , M and M 8 . 

In QCD, chiral symmetry is broken by the light quark masses. We will make 
the simplifying assumption, at the level of the sigma-model, to consider symmetry 
breaking to linear order in the quark masses but we accept all the terms which are 
soft (i.e. of dimensionality strictly smaller than four). Under these assumptions the 
most general symmetry breaking sector has, again, four parameters, 

C SB = 7o(Sx f ) + 7 i(S- 1 x)det(S) + 72 (££ t £x t > 

+ 7 3<££t><£ X t} + / i . c . . (3.22) 

Usually, only the first term is considered. This increased phenomenological flexibility 
will allow us to obtain better fits of the scalar nonet masses. While more general, this 
lagrangian will prove nevertheless to be reasonably constraining. We note that only 
the first two terms are renormalizable in the strict sense that the counterterms are 
exactly of the same form[|n]]. The last two terms generate counterterms which are of 
higher order in the scalar source x- Consistent with the assumption of not including 
further terms quadratic in the quark mass matrix at this level, we will compute the 
scalar meson masses and the vacuum expectation values at linear order as well. We 
have checked this approximation by comparing with exactly calculated masses. 

The connection with the representation of the scalar fields used in sec. 3.1 is 
performed by making a change of variables |3! 




£ = exp I i\j — -^r I uHu (3.23) 



3Fr 



o 



with 



« = e*p(*£^f), H=(v + ^=)I d + S . (3.24) 



This change of variable makes sense only if Fq ^ 0, i.e. when chiral symmetry is 
spontaneously broken. The kinetic energy part of the linear sigma-model lagrangian 
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gets transformed into 



-Jl^(H^) (3.25) 



with 



v^s = o^s + [r M , s\, r M = h^%u + ud^), = - ud^). (3.26) 

Replacing if by its expression above ( |3.24 ) one identifies v with the pion decay 
constant in the chiral limit Fq and one obtains a prediction for the two couplings Cd 

and dd, 

v 2 = y , c d = v, Cd = ^- (3.27) 
The symmetry breaking lagrangian, expressed in terms of the new fields reads 

Csb = ((loH + l2 H* + l3 (H) 2 H)( X+ cos ^ - i X - sin ^|^)> 

+ll detH(H-\ X+ cos J?^. + i X - sin £)) (3.28) 

In this sector, we obtain one relation among the four parameters 7* from the require- 
ment that the coefficient of the chiral lagrangian term (x+) be correctly normalized 
to Fq/4 and we can express two other parameters in terms of the meson-source 
couplings c m and c m . This gives the relations 

3 12 2 

70 = -V + r C ™ J^Cm + 73^ 

4 6 V3 

liv = -g(cm - v^c™) - 273W 2 

7 2 w 2 = — -« + -Cm + ~ 2 T3^ 2 • (3-29) 

4 6 V3 

In this sector, 73 as arbitrary parameters. 

In this model, we can now express the many parameters that appeared in the 
general scalar lagrangian discussed above. Firstly, for the terms coupling one scalar 
field to two scalar sources we have simply a' — V = d — d! = 0. For the trilinear 
couplings of resonances, one obtains the relations 

a = \=- ( Ml + -Ml 

b = -=- (Ml + 2M 2 - -M 2 

' ■ (3-30) 
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Finally, the parameters g m , h m , e m , f m which control linear symmetry breaking of the 
scalar meson masses (see eqs.( |3.12|) , ( |3.13| )) obey the following relations, 

-c m + c mJ ^)- l -8 l3 v 




V3/ m - - ^V3c m - c m ^ j-- + -g ri 



Ml - ZMl ( ^ MJ\ _ 3_ _ 1 



9M 2 « l V m m M|/ 16 



fc» = ~^-f ^-^-A-l^. (3.31) 
18M 2 w V M|y 32 24 y v ; 



Making use of the relations ( p3.11| ) we can express the three 0(p 6 ) LEC's Ci 9 , C 2 q, C 2 i 



discussed above in terms of the parameters of the linear sigma-model. We are also 
interested in the prediction for the OZI suppressed 0(p 6 ) LEC C\q which participates 
in the flavour variation of the order parameter F n as discussed above. After a small 
calculation, the following expression is obtained, 



F 2 C 16 = i(2c w + 10V3c 
54 




(3.32) 



3.3 Phenomenological applications 

Altogether, the sigma-model as considered here has six independent parameters 
(apart from Fq): three chiral limit masses Mq, Mg and M Vo and three couplings 
to the scalar source c m , c m and 73. The parameter M Vo , i.e. the mass of the rj' in the 
chiral limit can be estimated to be|Q, 

M Vo ~ 0.8 GeV . (3.33) 

Let us now discuss the determination of the remaining five parameters. A priori, 
one expects to be able to reproduce exactly the four independent masses (M ao , M K0 , 
M ao , Mf ) in the scalar nonet and, as an additional constraint a natural choice would 
be to enforce the super-convergence relation (|2.15|) which implies here c m = V3c rn . 



In practice, however, because of non-linearities, this choice strongly restricts the 
possible range of masses for the k . This is an indication that this model should be 
considered as a toy model rather than a really physically meaningful one. Instead, 
we did not try to impose the super- convergence relation and use a constraint from 
the pseudo-scalar sector, the ratio Fk/F^ which allows for a fairly broad range of 
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values for M KQ . Because of the specific parametrisation of the mass matrix, eq. ( 3.13|) 
its entries get evaluated as follows in terms of physical masses, 



M 



22 



4M K 2 o 



Ml 



M 
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M 2 o + M 2 fo - M 22 , 



(3.34) 



and, using the determinant, 

M 12 = ± 



-(M22-MI )(M 22 



(3.35) 



We observe that there are two possible sign assignments. Reality of M.\ 2 requires 
that the following inequality be satisfied 



< + 3 < < Ml < < +3M ^ , (3.36) 
4 — K ° — 4 

As shown in [S3], this inequality holds in a model independent way if symmetry 
breaking is assumed to be linear. From the expression of the mass matrix (|3.13|) it is 
then easy to see that the scalar singlet chiral mass M5 2 satisfies a quadratic equation 



M n 4 - BM 2 + CM. 



no 



(3.37) 



with 



B 

C 



1 



9 

-(2m 



mz 



1 



(2m 2 K + mi) 



' 1 + i< 



9 m 2 K 



1 M 2 — M 2 \ 
27 m 2 K — m 2 , J 

■K 1 ^ 



12 



3 



(3.38) 



Real solutions exist provided 



B 2 - ACM 2 > 



'/n 



(3.39) 



and it must further be checked that at least one solution is positive. This puts further 
constraints on the allowed values of the scalar meson masses in this model. Once Mq 
is determined, the other parameters are easily evaluated. For instance Mf is given 
by 



Ml 



1 



" 1 + ^ 
+ \{2m 2 K 



M z + 



2CM 



'/<> 



2 -^rn 2 K + mi) 



Ml 



rn 



K 



mi 



(3.40) 



Next, using as experimental input M 2 o — M„ and .M12 together with eqs.( |3.12|) , 
( p,13| ) and ( |3.31| )) gives two linear equations for the three quantities c m /v, c m /v and 
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73^- We will use as an additional constraint, the ratio Fx/F n , which determines 
c m /v from the relation, 

c m Mj F K -F T ^ 
v 2F n m\ - ml ' 

Finally, there remains to determine v = F / \[2 . To linear order in the quark masses, 
v is given by the following expression 



V2 



2 2 (2c m /v V3c m /v\ 4 2 f-c m /v VSc m /v' 
1 - . + — — T5— - -m K — = + 



m| m 2 ; 3 K V M l M o 



(3.42) 



This expression can be recovered in two different ways: one can either use CHPT 
together with the expressions ( |3.3| ) for L4 and L5 or write down the equations for the 
vacuum expectation values in the presence of quark masses and solve these equations 
to linear order in the quark masses. 

3.3.1 Light cr meson 

The existence of a very broad scalar resonance in inr scattering with M ~ T ~ 
0.5 — 0.6 GeV is by now well established since the work of Basdevant Froggatt and 
Petersen [1^1, in which the whole set of constraints from unitarity, analyticity, and 



crossing symmetry has been implemented (see ref. |35| for a recent survey and the 
particle data book|J for a complete list of references). What is unclear is whether 
this state should be interpreted as a physical light scalar resonance. Recently , Black 



et al.|20| have proposed arguments based on perturbative unitarity favouring the 
existence of a light a and also, of a light k meson. According to the inequality ( |3.36| ), 
if the sigma meson is light then one must have M K0 < M ao unless our assumption 
of linearity in the quark masses fails. Recently, Tornqvist |Hj (similar fits were also 
discussed in ref. [15[]) has attempted to accommodate in the linear sigma-model a 



scalar meson multiplet with a light sigma and a heavy k. Such fits fail to obtain 
the /o(980) at the correct mass and therefore would not correctly reproduce the nir 
phase-shifts around 1 GeV, even if one- loop corrections are included ||34||. 

Let us now assume that a light sigma meson exists, e.g. M UQ ~ 0.6 GeV, and 
discuss the consequences. We take otherwise M ao = 0.983 GeV, and Mf = 0.980 
GeV from experiment and M K0 = 0.9 GeV as proposed in ref. |PD" . From this input 



one easily calculates the parameters A and B in the quadratic equation ( |3.37| ), 

A = 0.024 GeV 2 B = 0.042 GeV 2 . (3.43) 

There results that no real solution to the equation for Mq exists unless M m < 0.06 
GeV, which appears as an absurdly small value. Reality and positivity of Mq (and 
M|) are necessary conditions for the existence of a stable minimum of the potential, 
with spontaneous chiral symmetry breaking. One must thus conclude that with such 
input scalar meson masses Np = 3 chiral symmetry is not spontaneously broken in 
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M. 


M 


M s 


Cm 


V3c m 


73 


V2v 


1.10 


1.17 


0.85 


0.023 


0.040 


0.54 


0.092 


1.20 


1.05 


0.93 


0.028 


0.027 


0.33 


0.094 


1.30 


0.87 


0.90 


0.023 


0.034 


1.56 


0.081 


1.33 


0.79 


0.85 


0.017 


0.037 


2.79 


0.068 


1.35 


0.72 


0.80 


0.011 


0.033 


4.82 


0.051 



Table 3: Determination of parameters in the case of a heavy scalar nonet M ao = 0.980, 
Mf = 1.5 for several values of M KQ . All entries are in GeV except 73 which is in GeV^ 1 . 



the linear sigma-model. The same conclusion holds if one picks up larger values for 
M KQ : using the set of parameters determined in ref.[|II| for M a = 0.4 and M a = 0.6 
GeV one verifies that the equation for the vacuum expectation value v eq.( |3.18 ) has 



no real solution v 7^ 0. A better situation prevails if we accept smaller values for M KQ . 
If we take, for instance, M KQ = 0.750 GeV, which is the smallest value allowed by the 
inequality ( |3.36| ), then, eq.( |3.37| ) has real solutions provided M m < 0.62 GeV, which 



is still reasonable. A problem remains, however, that being close to a situation where 
chiral symmetry is unbroken for Np = 3, while it is in general broken for Np = 2, 
the expansion in the strange quark mass is likely to be non-converging. Indeed, if we 
try to compute F with, say, M m = 0.50 GeV, using our linear expansion formula 
( |3.42 ), we find Fq ~ 28 MeV, which is very small compared to F n , suggesting that 



the expansion in m s is not perturbative in this case. 
3.3.2 Heavy a meson 

Let us now consider the scenario of a heavy sigma meson, i.e. we identify the sigma 
meson with the well established resonance /o(980) and we identify the heavier 1 = 
member of the nonet with the /o(1500) (see pl[ ). With these assignments, we find 



that a solution of eq. (|3.37|) with real Mq exists in the range 

M K0 < 1.37 GeV , (3.44) 

taking M m = 0.8 GeV and the negative sign for M.12 (see ( |3.35| )). This means that 
we cannot exactly reproduce the experimental mass of the Xq(1430) but we can get 
reasonably close. We will eventually allow M Ko to vary somewhat away from the 
experimental result. Numerical values of the parameters M , M 8 , c m , c m , 73 as well 
as v determined for several values of M K0 are collected in table (3). We note that 
the situation with a minimal symmetry breaking lagrangian (i.e. 71 = 72 = 73 = 
) is close to the case when M Ko = 1.20 in the table. In this case, one has c m = 
\[?>c m = v/2. The value of F = \/2v remains reasonable close to F n except when 
M KQ gets very near to its upper bound. The coupling c m comes out smaller than 
that determined in ref.[(| (see ( |3.5| )) while the coupling Cd = v comes out larger. The 
coupling q controls the decay width a — > nr], which tends to be too large in this 
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M K0 


10 3 L 5 


10 3 L 7 


10 3 L 8 


io 5 c 19 


io 5 c 20 


1.10 


2.06 


-0.22 


0.36 


-0.17 


-0.05 


1.20 


2.15 


-0.47 


0.45 


-0.34 


0.07 


1.30 


1.61 


-0.62 


0.32 


-0.45 


0.12 


1.33 


1.13 


-0.51 


0.20 


-0.38 


0.08 


1.35 


0.64 


-0.31 


0.10 


-0.24 


0.04 



Table 4: Numerical values of some 0(p 4 ) and 0(p 6 ) LEC's from the sigma-model with a 
heavy sigma. 



model (accepting that the experimental value is T ~ 60 MeV, which is subject to 
some debate p), the result in this respect improves as M K0 gets larger. 

Once these couplings are known it is easy to determine the LEC's with the 
formulas given above. Numerical values for the LEC's dominated by the scalar 
mesons and not suppressed in the large N c limit are shown in table (4). We have 
also included L 7 , which is assumed to be saturated by the r] , computed from the 
expression [||, 

vo 

where the coupling d m of the r/o field to the pseudo-scalar current has the following 
expression in the present version of the linear sigma-model, 

~ dm = ~276 v 1 + ~ ^~ Cm) + Ul3V ) ' (3 ' 46) 

The fact that L5 does not remain constant when M Ko increases is a reflection of the 
effect of non-linearities in m s . The values of the LEC's may be compared with the 
ones determined at 0(p 4 ) from experimental data|I[ 

lO 3 ^^) = 2.2 ±0.5, 10 3 L 7 = -0.4 ±0.15, 10 3 L r 8 (m v ) = 1.1 ± 0.3 . (3.47) 

While the prediction for L 7 is acceptable, that for L 8 is too small. This is related 
to the small size of the meson to source coupling c m predicted by the model. Still, 
the order of magnitude is correct except, perhaps, for the very last line in the table. 
The LEC's were recently re-determined in an 0(p 6 ) analysis of the experimental 



data|36fl and this brings the values somewhat down in magnitude: 10 /^(m^) = 
1.45±0.12, 10 3 L 7 = -0.31±0.15, lO 3 ^^) = 0.68 ±0.18. A resonance-saturated 
estimate for the 0(p 6 ) coupling C19 was provided in ref.0]. The value that we 
obtain is significantly smaller, by a factor 3 — 4. One difference in the estimates is 
that we include the effect of cubic interaction terms, but the numerical influence of 
these terms turns out to be unimportant. The main reason for our smaller result is, 
again, the fact that the coupling c m is smaller here. 
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M, 


10 3 L A 


10 3 L 6 


10 5 C 16 


10 5 (C 20 + 3C 21 ) 


1.10 


-0.06 


0.07 


0.19 


0.02 


1.20 


-0.16 


-0.04 


0.08 


0.05 


1.30 


0.33 


0.15 


0.21 


0.08 


1.33 


0.56 


0.28 


0.36 


0.15 


1.35 


0.55 


0.31 


0.41 


0.27 



Table 5: Numerical values of OZI suppressed 0(p A ) and 0(p 6 ) LEC's from the sigma- 
model with a heavy sigma. 

Let us now return to the OZI suppressed coupling constants L 4 and L 6 and their 
0{jP) counterparts C\q and C20 + 3C21. The predictions from the linear sigma-model 
are collected in table (5). The values and signs of L 4 , L 6 are seen to depend very 
much on the mass of Kq. They both approximately vanish when M K0 ~ 1.24 GeV. 
This is the point of minimal OZI violation. As one pushes the k,q to higher masses, 
as required by experiment, then the pattern of OZI violation is not unlike the one 
found on the basis of sum rules. Both L 4 and L 6 become positive, and the orders of 
magnitudes seem to be in agreement with the discussion based on sum rules. There 
is even a relatively good agreement as far as L 4 and Ci 6 are concerned. 

4. Conclusions 

In this paper, we have pursued an investigation on the evolution of order parameters 
of the QCD chiral vacuum as one increases the number of massless flavours, as can 
be inferred from experimental data. An information on this evolution is contained in 
the 0(p 4 ) coupling constants L 4 and L 6 which control how F n and (uu) respectively 
vary from the SU(2) chiral limit with m s ~ 200 MeV to the £77(3) chiral limit. 
Estimates for the couplings L 4 and L 6 are based on the scalar form factors of the 
pion and the Kaon which are reconstructed from experimental S-wave scattering 
data modulo some assumptions. We have investigated the 0{m s ) corrections to our 
previous results. These have two origins: a) the normalization of the form factors at 
the origin must include the 0(p 4 ) contributions and b) the chiral expansions (for F n 
and II 6 (0) ) must be pursued up to order p 6 . 

These corrections seem not to affect in a significant way the results obtained 
previously. In particular there is no sign that these corrections go in the sense 
of decreasing the values of L 4 and L§ to make them compatible with the naive 
large N c expectation. This conclusion, however, only holds provided one makes the 
assumption that the chiral expansion is reasonable, in other terms, that the 0{p & ) 
contribution (say in the expansion of II 6 (0)) is not larger than the 0{p 4 ) contribution. 
Modulo this assumption, one also obtains from this analysis an estimate of the 0(p e ) 
coupling constant C±q and the combination C20 + 3C21. Influence of the values of L 4 , 
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Lq on the convergence rate of the chiral expansion for various observables was studied 
very recently in [0. Because of such assumptions which enter into the calculation, 
the error bars on our coupling constant estimates must be considered as educated 
guesses. Our results are compatible with a tendency towards chiral restoration upon 
increasing Np from Np = 2 to Np = 3 which is significant, particularly so for the 
quark condensate (in agreement with the arguments of 0), for which the decrease 
by a factor of approximately two found in ref.0 is confirmed. This behaviour of 
the quark condensate is in qualitative agreement with the one obtained from the 
instanton liquid model |38]. It would be also interesting to compare this result with 
unquenched lattice simulations with Np = 3. Unfortunately, such simulations are 
not available yet, but will exist in the near future. 

We also discussed some predictions from the linear sigma-model, considered as 
a simple toy model. One application is to gauge the possible importance of self- 
couplings of resonances (which are present in the model) in resonance saturation 
estimates of the 0(p 6 ) LEC's. Furthermore, the model provides a connection between 
simply the spectrum of the scalars and the evolution of chiral order parameters with 
Np. In order to improve the ability of the model to reproduce experimental masses, 
the symmetry breaking sector was made more general than usually done, but still 
linear in the quark mass matrix. We have considered two possibilities for the light 
scalar nonet: 

a) The nonet contains a light a and a light k. In this case, we find that chiral 
symmetry tends to be restored already at Np = 3 or, at least, for a value of Np 
somewhat too small, making the chiral expansion in m s unreliable. 

b) The nonet is composed of the resonances a (980), / (980), i^g(1430) and / (1500) 
(which are well established). In this case, we find that we need only a small amount of 
nonlinear effects to exactly reproduce the -K"q(1430) mass, and SU(3) chiral symmetry 
is realised in the Goldstone mode. Concerning L4 and Lq, we find that they become 
different from zero as one increases the mass of the Kq(1430) towards its experimental 
value. The signs and orders of magnitude, then, are compatible with the preceding 
analysis. The same is true of the related 0(p 6 ) constants. 
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A. The correlation function II 6 (s) to 0(p 6 ) 

At chiral order p 4 , a simple calculation gives the correlation function Il 6 (s) 



n 6 (s) = 2J K (s) + p ri (s) + 6AL r 6 



1 



16tt 2 



4 

2-^^ + -L v 
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where Jp(s) is the one- loop function defined as to vanish for s = 

1 



Jp(s) 



16tt 2 



Op — 1 
crp log — r + 2 ) , a P 



4Mp 



(A.I) 



(A.2) 



(7p + l 

In all the formulas Mp stands not for the physical pseudo-scalar meson masses but 
for their lowest order chiral expansion, 

M* = 2mB , M* K = (m + m s )B , M;= (2m + ^ . (A.3) 
The following simplified notation for logarithms was introduced 



log 



Ml 
V 2 



(A.4) 



At order p 6 , one must compute the one- loop and two-loops diagrams shown in fig.l 
and add the tree-level contributions from the 0(p G ) chiral lagrangian. The calculation 
is lengthy but not excessively difficult because of the absence of sunset-type diagrams 
in the present case. The result may be written as a sum of six terms: 

An 6 ( s ) = An^) + An^( s ) + Arr 6 esc ( s ) + a-Jj + b + c . (A.5) 

The first term encodes the 0(p 4 ) corrections to the pseudo-scalar meson masses in 
the one-loop functions, 
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(A.6) 



These mass corrections are given by the following expressions! 
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In these expressions, 5*6, S 7 and S 8 denote the following combinations of 0{p 4 ) cou- 
pling constants, 
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The second corrective term An|(s) is generated by 0(p A ) corrections to the 
scalar form factors and has the following form, 
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The derivatives of AM|- and AM^ are easily computed from the formulas above and 
for AMj one has, 
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The third corrective term AITg esc (s) is the rescattering contribution, 
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The last corrective term is a polynomial, linear in s. This term picks up contributions 
from four coupling constants of the 0(p 6 ) chiral lag rangian, C39, C20' 21 <md Cg4 in 
the notation of [|TU |. The part proportional to s is as follows 
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The constant terms, finally are 
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(A.15) 



The terms proportional to the mass corrections AM|- and AM 2 (see eg. ( |A.6| ) and 
( |A.15| )) may eventually be absorbed into the 0(p 4 ) expression for Il6(s), which 
amounts to replace the lowest order expression for the pseudo-scalar masses there 
by their expression up to 0(p 4 ). Formula (|2.23|) in the text is obtained from B + C, 
setting m = 0. 
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Figure 1: Set of one and two-loops Feynman graphs contributing at chiral order p 6 to the 
correlation function Il6(s): crosses denote vertices from the 0(p 2 ) chiral lagrangian, boxes 
denote vertices from the 0(p 4 ) chiral lagrangian. Some diagrams obtained by exchanging 
external lines have not been drawn. 
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